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ABSTRACT 

J 

For  unsteady,  irrotational  flow  governed  by  Laplace fs 
equation,  velocity  potential  and  stream  function  solutions  are 
presented  with  particular  consideration  being  given  to  the  boundary 
conditions.  The  analyses  are  applicable  to  multiple  bodies  each 
moving  in  an  arbitrary  direction  with  varying  velocity.  In  the  final 
section  of  the  report,  the  problem  of  calculating  the  entrained  mass 
for  a  body  of  arbitrary  shape  is  considered. 


1. 


1.  INTRODUCTION 

In  previous  reports  [1,2]  the  application  of  the  finite 
element  technique  to  steady,  irrotational,  incompressible  flow  fields 
was  considered.  The  special  boundary  conditions  encountered  in 
potential  flow  problems,  such  as  the  Kutta  condition  on  an  aerofoil 
were  dealt  with  in  some  detail. 

This  present  report  extends  this  consideration  to  unsteady 
flow  fields,  concentrating  particularly  on  potential  flow  problems. 

The  scope  is  restricted  to  two  dimensions.  The  extension  to  three 
dimensions  is  under  investigation. 

2.  UNSTEADY  FLOW  ANALYSIS 

In  an  unsteady  field,  the  boundary  conditions  vary  with  time. 
Since  potential  fields  are  uniquely  specified  by  the  boundary  conditions, 
there  is  a  corresponding  variation  with  time  in  the  field.  The  potential 
at  a  point  in  the  domain  changes  as  signals  reach  that  position,  from  the 
boundaries.  The  velocity  of  these  signals  is  the  acoustic  velocity, 
which  in  an  ideal  fluid,  is  infinite  in  magnitude.  Thus  for  potential 
flow,  changes  in  the  boundary  conditions  affect  immediately  all  points 
in  the  field.  The  stream  function  and  velocity  potential  distributions 
are  thus  always  in  accord  with  the  boundary  conditions  at  any  given 
instant.  This  does  not  mean,  however,  that  steady  state  conditions 
exist  at  every  instant  of  time.  The  ideal  fluid  has  inertia,  and  there 
is  thus  a  corresponding  acceleration  term  which  must  be  included  in  the 
equations  of  motion.  The  solution  of  these  equations  will  thus  be 
different  from  that  obtained  using  the  steady-state  equations,  even  for 


identical  boundary  conditions. 


As  in  steady  flow  problems,  it  is  desirable  to  formulate  the 
field  equation  in  terms  of  a  single  scalar  variable.  This  is  simply 
done  by  means  of  the  velocity  potential  or  the  stream  function.  There 
are  some  difficulties  in  the  boundary  conditions  that  arise  if  the 
stream  function  is  used,  and  these  are  discussed  in  Section  3.2.  For 
constant-density  flows  with  vorticity,  an  acceleration  potential  can 
be  derived  which  satisfies  the  Poisson  field  equation.  This  approach 
is  particularly  valuable  since  it  can  be  applied  to  rotational  flows. 
This  report,  however,  deals  only  with  the  case  of  the  irrotational 
field. 


2.1  The  Velocity  Potential 

For  unsteady , irrotational  flow,  a  velocity  potential  <p  can 
be  defined  from  the  consideration  of  zero  vorticity  as  follows. 

The  vorticity  \  is  defined  by 

l  -  V  x  q  ,  (2.1) 

where  q  is  the  velocity  vector,  which  in  terms  of  its  x,y,  components, 
u,v,  respectively,  is  given  by 

q  =  ui  +  vj  ,  (2.2) 

and  where  V  is  the  Laplace  operator.  In  Eq.  (2.2),  i  and  J  are  the 
unit  base  vectors  in  the  x  and  y  direction,  respectively. 

The  condition  of  zero  vorticity,  i.e.. 


3. 


becomes  from  Eq.  (2,1) 

V  x  q  =  0  .  (2.4) 

Equation  (2.4)  may  be  satisfied  by  defining  the  velocity  vector,  q,  in 
terms  of  a  scalar  $  as* 

q  =  -  V<j>  .  (2.5) 

The  scalar  <J>  in  Eq.  (2.5)  is  known  as  the  velocity  potential.  From 
Eq.  (2.5),  the  velocity  components  u,v,  respectively,  become 


u 


3j) 

3x 


(2.6a) 


v 


2& 

3y 


(2.6b) 


If,  in  addition,  the  fluid  is  assumed  to  be  incompressible, 
then  it  can  be  shown  [3,4]  that 

V  *q  =  0  .  (2.7) 

Substituting  Eq.  (2.5)  into  Eq,  (2.7)  results  in 

-V.(Vcf ))  =  0  ,  (2.8a) 

which  further  reduces  to 

v2<t>  =  0  .  (2.8b) 

The  combination  of  the  irrotationality  and  continuity  con¬ 
ditions  therefore  yields  Laplace’s  equation  in  $ .  Although  the 


*  The  convention  that  q  *  -  V<f>  is  adopted  in  this  report.  However, 
some  authors  prefer  the  convention  q  *  +  V<p. 
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independent  variable  t  may  appear  in  the  expressions  for  the  velocity 
vector  q ,  and  consequently  in  the  velocity  potential  <f>,  the  above 
analysis  shows  that  the  field  equation  for  the  unsteady,  incompressible, 
irrotational  flow  is  the  same  as  that  for  the  steady  case.  The 
boundary  conditions  are,  however,  the  appropriate  unsteady  ones,  and 
are  dealt  with  in  Section  3.1. 

The  pressure  at  any  point  in  the  fluid,  p,  may  be  obtained 
from  the  generalized  Bernoulli! ’s  equation  [3,4] 

H  +  !  (V<t>)2  +  ^  *  f(t)  ,  (2.9) 

where  p  is  the  density  of  the  fluid,  and  f(t)  is  a  space  constant.  It 
is  noted  that  for  the  unsteady  case  both  space  and  time  derivatives  of 
the  velocity  potential  <f>  are  required.  The  space  constant  f(t)  in 
Eq.  (2.9)  allows  for  a  change  in  the  absolute  pressure  by  an  arbitrary 
variation  of  the  external  pressure  on  the  system.  If  the  external 
pressure  Is  considered  fixed,  then  f(t)  reduces  to  a  simple  constant, 
both  in  space  and  time. 

Although  the  velocity  potential  <f>  is  of  importance  in  the 
analysis  of  both  steady  and  unsteady  flows,  it  is  rather  restrictive 
in  that  it  can  only  be  used  if  the  flow  is  irrotational. 

2 . 2  The  Stream  Function 

Laplace’s  equation  similarly  applies  to  the  stream  function 
^  in  unsteady,  irrotational  flow,  if  the  fluid  is  incompressible.  This 


may  be  shown  as  follows. 
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The  condition  of  incompressibility,  i.e,  Eq.  (2.7), 

V.q  =  0  ,  (2.7) 

may  be  satisfied  by  choosing  the  velocity  vector  q  as 


From  Eq. 


q  = 

(2.10) 


|*I  +  |*j 

3y  3x  J 


,  the  velocity  components  u,v,  respectively,  become 


(2.10) 


u 


d\p 

3y 


(2 


v 


(2.11b) 


Using  the  condition  of  irrotationality ,  Eq .  (2.4),  there  is 

obtained 


lli  +  iii 

9x2  T  9y 2 


0 


9 


which  is  Laplace’s  equation  in  0,  i.e., 
V2^  =  0 


(2.12a) 


(2.12b) 


The  statement  made  regarding  the  field  equation  for  the  velocity 
potential  <p  in  Section  2.1  similarly  applies  to  the  stream  function  ^ , 
irregardless  of  whether  the  independent  variable  t  appears  explicitly  or 
not  in  the  velocity  vector  q  and  (consequently)  in  the  stream  function 
\p. 

The  finite  element  technique  may  be  applied  to  solve  for  ij;  in 
the  manner  described  in  Section  4.2.  The  boundary  conditions  will 
generally  be  more  complex  than  those  for  the  velocity  potential.  These 


6. 


conditions  are  considered  in  detail  in  Section  3.2. 


2 . 3  The  Cauchy-Riemann  Conditions 

It  was  shown  in  Sections  2.1  and  2.2  that  for  unsteady,  in¬ 
compressible,  irrotational  flow,  both  a  velocity  potential  $> ,  and  a 
stream  function  \p  exist,  and  furthermore  that  both  these  scalar  fields 
<p  and  ^  satisfy  Laplace’s  equation. 

It  is  instructive  to  note  that  although  the  unsteady  flow  is 
being  considered,  i.e., 

<p  =  <f>(x,y,t)  >  (2.13a) 

ip  *  ^(x,y,t)  ,  (2.13b) 

the  Cauchy-Riemann  conditions  still  apply.  This  may  be  shown  as  follows, 
where  throughout,  the  explicit  dependence  of  the  dependent  variables 
q,  u,  v,  <j>,  and  \p  on  the  independent  variables  x,  y,  and  t,  is  included. 
From  Eq.  (2.2)  there  is  obtained 

q(x,y , t)  =  u(x,y,t)i  +  v(x,y,t)j  .  (2.14) 

Equations  (2.6a)  and  (2.6b)  similarly  reduce  to 


u(x,y,t) 


34>(x,y  ,t) 
9x 


(2.15a) 


and 


v(x,y,t) 


a$(x»y»t) 

dy 


(2.15b) 


Furthermore,  from  Eqs.  (2.11a)  and  (2.11b),  it  follows  that 


u(x,y,t) 


v(x,y,t) 


(2.16a) 


(2.16b) 


Comparing  the  set  of  equations  (2.15)  with  the  set  (2.16)  results  in 


(2.17a) 


(2.17b) 


which  shows  that  although  <p  and  are  dependent  on  the  variable  t,  they 
still  satisfy  the  Cauchy-Riemann  conditions.  The  above  analysis  can  be 
extended  to  show  that,  in  addition  to  the  set  of  Eqs.  (2.17),  the 
following  relations  hold. 


(2.18a) 


3<fr(x,y,t) 

3s 


(2.18b) 


where  -  and  -  are  the  velocities  in  the  directions  of  increasing 
Bn  Bs 

n  and  s,  respectively.  The  vectors  n  and  s  are  mutually  orthogonal  and 
may,  for  example,  be  chosen  as  the  unit  outward  normal  to  the  curve  S, 
and  as  the  unit  vector  tangential  to  S,  respectively,  see  Figure  2.1. 


< 


Velocity: 
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Fig.  2.1  Velocity  components  in  the  directions  n  and  s. 

2.4  The  Acceleration  Potential 

A  potential  function  P  can  be  defined  for  acceleration,  where 
the  flow  is  such  that  the  vorticity  is  preserved  [3].  Its  special 
importance  is  that  it  may  be  defined  for  both  rotational  and 
irrotational  flows,  whereas  the  velocity  potential  can  be  defined  for 
the  latter  case  only,  as  had  already  been  mentioned  in  the  preceding 
Section  2.1.  The  relevant  field  equation  in  terms  of  P  will  not  be 
derived  here,  but  it  can  be  shown  to  satisfy  Poisson's  equation.  The 
analysis  in  terms  of  the  acceleration  potential  P  is  especially  useful 
in  the  analysis  of  discontinuous  flews  since,  although  the  velocity  is 
discontinuous  across  a  vortex  sheet,  and  the  velocity  potential  <p  may 
also  be  discontinuous  there,  the  pressure  and  hence  the  acceleration 
potential  P  are  continuous  across  such  a  sheet. 

The  authors  are  presently  investigating  the  application  of 
the  finite  element  technique  to  such  discontinuous  flows. 
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3-  BOUNDARY  CONDITIONS 

For  a  stationary  impervious  boundary,  the  boundary  condition 
is  of  the  homogeneous  Neumann  type  [  5,  6  ],  but  for  a  moving  boundary 
(whether  that  of  a  rigid  surface,  a  deforming  surface,  or  a  two-phase 
surface)  the  kinematic  boundary  condition  will  be  of  the  nonhomogeneous 
Neumann  type  as  is  shown  in  the  subsequent  sections.  In  the  following, 
only  irrotational  flow  is  considered. 

3.1  The  Velocity  Potential 

In  considering  the  boundary  conditions  for  the  velocity 
potential,  it  is  best  to  consider  a  particular  example.  The  same 
example  will  be  used  again  when  dealing  with  the  boundary  conditions 
for  the  stream  function  \p  in  Section  3.2. 

The  example  of  interest  is  that  as  shown  in  Figure  3.1,  namely 
two  rigid  objects  in  motion^  in  an  otherwise  uniform  stream. 


t  In  this  and  the  following  analysis,  the  velocities  of  the  objects 
may  depend  explicitly  on  the  time  t. 
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Fig.  3.1  Two  moving  objects  in  an  otherwise  uniform  stream. 


Let  the  two  objects  be  denoted  by  and  as  shown.  Furthermore  let 
the  boundaries  of  A^,  A^  be  given  by  the  (closed)  curves  S^, 
respectively,  and  let  their  corresponding  velocities  be  given  by  u^  and 
u^*  It  is  assumed  in  the  analysis  presented  throughout  this  report 
that  the  objects  are  impervious,  consequently  there  can  be  no  flow 
across  the  boundaries  and 

For  purposes  of  illustration,  the  following  analysis  only 
considers  one  object,  which  can  then  be  generalized  to  two  or  more 
objects  moving  in  the  same  flow  field.  Consider  therefore  the  case 
illustrated  in  Figure  3.2,  where  the  object,  denoted  by  A^  and  bounded 


by  S^,  is  shown  at  its  position  at  time  t,  i.e.,  with  respect  to  the 


fixed  coordinate  frame  Oxy. 


uniform 

stream 


uniform 

stream 


Fig.  3.2  Object  moving  in  an  otherwise  uniform  stream. 


Let  the  outer  boundary  of  the  control  volume  being  considered  be  denoted 
by  Sq,  which  is  assumed  to  be  sufficiently  far  removed  from  A^,  and  hence 
from  such  that  any  movement  of  the  object  A^  does  not  affect  the 
uniform  stream  exterior  to  Sq,  Acceleration  of  the  object  A^ 
is  allowable  in  this  analysis,  and  hence  the  velocity  of  A^  is  indicated 
as  u^(t),  see  Figure  3.2. 

The  finite  element  analysis  presented  in  this  report 
considers  the  solution  at  fixed  instants  of  time.  Suppose  the  solution 
is  sought  at  time  t,  i.e.,  when  the  object  is,  with  respect  to  the 
chosen  axes  Oxy,  at  its  location  as  shown  in  Figure  3.2.  Let  the  unit 
outward  normal  to  be  denoted  by  then  from  the  condition  that  no 
flow  occur  across  this  boundary  (since  S,  is  impervious) ,  it  is  clear 
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that  the  normal  velocity  of  the  boundary  of  the  object  must  be  the  same 
as  the  normal  component  of  the  velocity  of  the  fluid  which  is  adjacent 
to  this  boundary.  The  velocity  of  the  object  in  the  direction  n^, 
which  will  be  denoted  by  u^n,  is  given  by 

**  ■  •  (3.1) 

If  the  velocity  of  the  fluid  in  this  same  direction  and  along  the  same 

boundary  S-  is  denoted  by  qcn,  then  it  follows  from  q  n  -  u^  that 
l  sx  S1  1 

qs  &  u^.n^  •  (3.2) 


From  Eqs.  (2.18a)  and  (2.18b),  it  is  known  that  q  n  may  also  be  written 

S1 

as 


on 


(3.3) 


Substitution  of  Eq.  (3.2)  into  Eq.  (3.3)  results  in 

!£  ■  •  «i*®i  °n  si  • 

It  is  important  to  note  that  the  boundary  condition  given  by  Eq.  (3. A) 

is  of  the  nonhomogeneous  Neumann  type.  This  is  in  contrast  to  the 

homogeneous  boundary  condition  *  0  on  S,,  which  occurs  when  A„  is 

dn  1  1 

stationary  [  6 ] . 

Since  the  flow  exterior  to  Sq  is  assumed  to  be  uniform  at  all 
times,  it  follows  that  on  S^  the  velocity  potential  <J>  Is  known,  at  least 
apart  from  an  arbitrary  constant  of  integration,  consequently  on  S^  a 
Dirichlet  condition  exists  for  all  times,  as  is  shown  below. 
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Since  the  flow  exterior  to,  as  well  as  on  Sq,  is  uniform,  it 


follows  that 


u(x,y,t)  =  U 


v(x,y,t)  -  V 


0  ’ 


(3.5) 


where  U,  V  are  the  magnitudes  of  the  uniform  velocity  in  the  x,y 
directions,  respectively.  From  Eqs.  (2.15a)  and  (2.15b)  it  follows  that 


(3.6a) 


li  -  .  v  (  °n  S° 

9y  ) 

and  furthermore,  the  condition  that  the  flow  is  uniform  gives 


with  solution 


♦  =  -  Ux  -  Vy  +  C  on  S. 


(3.6b) 


(3.7) 


where  C  is  an  arbitrary  constant  of  integration  which  may,  without  loss 
of  generality^*  be  chosen  to  be  zero.  Consequently  the  Dirichlet  con¬ 
dition  for  (p  on  Sq  reduces  to 


4>  =  -  Ux  -  Vy  on  S 


(3.8) 


The  boundary  condition  for  (p  on  could  have  been  derived 
equally  as  well  as  a  nonhomogeneous  Neumann  boundary  condition  as 


t  The  choice  C  ■  0  also  corresponds  to  a  change  in  position  of  the 
coordinate  frame  Oxy,  which  is  arbitrary. 


14. 


follows*  Noting  from  Eq.  (3.6b)  that  =  0,  it  follows  that 


ii  *  ii  i*  +  It  iz  on  s 

3n  3x  3n  3y  3n  0 


(3.9) 


(3.10a) 

(3.10b) 


Fig.  3.3  The  normal  and  tangential  vectors  n^  and  s^ 
on  the  boundary  S^. 


where  (n^)^  and  (n^)g  are  t*ie  X  an<*  ^  components  of  the  unit  outward 
normal  to  SA,  n  ,  Eq.  (3.9)  reduces  to 


15. 


By  using  the  set  of  equations  (3.6a),  Gq.  (3.11)  may  be 
reduced  further  to 


3£ 

3n 


-  U<\>0 


v<y0  on  so 


(3.12) 


It  is  always  possible  to  choose  the  coordinate  frame  Oxy  in 
such  a  way  that  either  U  or  V  is  zero,  and  furthermore  the  boundary 
Sq  may  be  chosen  as  a  rectangle  so  that  in  addition  to  U  or  V  being 
zero,  the  components  (n^)g  and  (n^^  ^ecome  simPly  either  +1  or  -1, 
which  simplifies  the  calculations.  In  passing  it  is  noted  that  for  an 
arbitrary  shaped  curve  S^,  the  corresponding  boundary  condition  to 
Eq.  (3.12)  can  be  written  as 


3$  —  — 

t1,  ■  -  q  on  S. 

Bn  ^uniform  0  0 


(3.13) 


where  <lun^form  *s  t*ie  velocity  vector  of  the  uniform  flow  exterior  to 


Sq ,  i.e. , 


^uniform  "  Ul  + 


(3.14) 


It  was  shown  previously  that  on  the  non-homogeneous  Neumann 
condition,  Eq.  (3.4),  must  be  satisfied.  If  the  velocity  of  the  object 
is  prescribed,  then  its  location,  as  well  as  its  unit  outward  normal 
n^f  is  prescribed  at  all  times  with  respect  to  Oxy.  Consequently  at 
every  instant  of  time  Eq.  (3.4)  is  determined,  and  hence  also  the 
boundary  condition  on  . 

In  summary  then,  the  flow  field,  in  terms  of  the  velocity 
potential  <f>,  may  be  obtained  for  any  time  t  as  the  solution  (see 


Section  2),  to 


Vz4> 


0  in  D 


(2.8b) 


subject  to  the  Dirichlet  condition 


4>  =  -  Ux  -  Vy  on  Sr 


(3.8) 


or  the  nonhomogeneous  Neumann  condition 


3  <j> 

3n 


-  q  .c  •  n.  on  S. 
uniform  0  0 


(3.13) 


and  the  nonhomogeneous  Neumann  condition 


3n 


-  ui  .  ni  on  8 


(3.4) 


It  is  noted  that  at  any  given  time  is  prescribed  and  the 


location  and  orientation  of  is  fixed  and  hence  n^  is  prescribed, 


consequently  the  boundary  condition  on  S^,  see  Eq.  (3.4),  is  completely 


prescribed.  For  additional  bodies,  Eq.  (3.4)  is  applicable  also. 


3.2  The  Stream  Function 

Consider  again  the  example  used  in  Section  3.1,  namely  the 
unsteady,  irrotational  flow  field  due  to  a  uniform  stream  of  an 
incompressible  fluid  obstructed  by  moving  objects,  as  shown  in  Figure 
3.2.  Arguments  similar  to  those  used  in  that  section  indicate  that  4^ 
is  completely  prescribed  on  Sq  at  least  up  to  an  arbitrary  constant  of 
integration,  which,  as  was  indicated  in  Section  3.1,  may  be  chosen  to 
be  zero.  This  procedure  is  illustrated  as  follows.  The  boundary 
condition  for  \p  on  Sn  is  obtained  from  the  conditions,  see  Eq.  (3.5), 


17. 


u(x,y,t)  -  U 


v(x,y ,t)  =  V 


on  S, 


(3.5) 


which  in  view  of  Eqs.  (2.16a)  and  (2.16b)  result  in 


* 


Vx  -  Uy  +  C  on  Sr 


(3.15) 


In  Eq.  (3.15)  the  arbitrary  constant  of  integration  is  denoted  by  C, 
which  may  be  chosen  to  be  equal  to  zero,  thus  reducing  Eq.  (3.15)  to 


ip  *  Vx  -  Uy  on  S 


(3.16) 


In  a  similar  fashion  to  that  discussed  for  the  velocity  potential  <J>, 
the  boundary  condition  for  \p  on  may  be  written  as  a  nonhomogeneous 
Neumann  condition.  This  may  be  shown  to  be 


H  ■  v(Vo  -  D(Vo  °"so 


(3.17a) 


or 


It 

3n 


^uniform  x  n0^ 


k  on  S 


0 


(3.17b) 


where  k  is  the  unit  base  vector  in  the  z  direction,  perpendicular  to 
both  i  and  J,  see  Figure  3. A. 


z 

Fig.  3.4  The  base  vectors  i,  j,  and  k. 

The  other  boundary  condition  to  be  investigated  is  that  on  the 

boundary  of  the  object  A^.  It  is  worth  reiterating  at  this  point 

* 

what  happens  on  this  boundary  in  the  steady  case.  In  the  steady  case, 
the  stream  function  \p  *  i|;(x,y)  may  simply  be  written  as 

ip  *  i|/ (s)  on  ,  (3.18) 

where  s  is  measured  along  from  some  arbitrary  datum,  i.e,,  from  some 

arbitrarily  chosen  point  P,  where  s  **  0.  The  condition  that  no  flow 
occur  across  the  impervious  boundary  dictates  that  the  velocity 
normal  to  must  be  zero.  Consequently,  from  Eq.  (2.18a)  there  obtains 

-  0  on  S,  ,  (3.19) 

O  S  1 

and  hence 


when  the  object  is  stationary. 
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^  *  constant  on  .  (3.20) 

It  needs  pointing  out,  that  although  ip  is  constant  on  in  the  steady 
case,  this  constant  is  not  known  in  value.  However,  it  may  be  obtained 
by  a  further  consideration  of  the  boundary  conditions  on  as  is  outlined 
in  the  following  sections,  as  well  as  in  the  literature  [2,7]. 

Extending  the  above  argument  to  unsteady  flow,  it  is  clear 
that  at  any  given  instant  of  time  t  the  stream  function  ip  must  still  be 
a  function  of  s,  i.e., 

ip  *  \p(s)  on  S1  ,  (3.21) 

at  that  particular  instant  of  time.  In  general,  it  would  be 
more  appropriate  to  make  the  statement  that 

ip  =  iKs.t)  on  S1  ,  (3.22) 

since  as  time  varies,  the  location  of  is  changing  accordingly.  It 
was  pointed  out  in  Section  3.1,  however,  that  in  the  present  approach 
the  finite  element  solution  is  to  be  obtained  at  a  fixed  instant  of 
time,  hence  the  formulation  (3.21)  will  suffice. 

Since  the  object  moves  with  time,  the  fluid  adjacent  to 
S^  must  have  a  velocity  associated  with  it  and  consequently  the 
expression  (3.19)  no  longer  holds.  From  Eq.  (3.21),  there  obtains 

di/>(s)  *  ^  ds  on  S^  .  (3.23) 

Integrating  Eq.  (3.23)  yields 
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<Ks)  *  i>(0)  + 


fids 

ds 


on  S, 


(3.24) 


Using  the  relation  (2.18a),  i.e.. 


3<{>  _ 

9n  3s 


(2.18a) 


and  noting  that  4^  =  for  the  case  being  considered,  Eq,  (3.24)  may 
ds  dS 


be  written  as 


Ms)  -  ’MO)  + 


f^-  ds  on  S 
3n  1 


(3.25) 


Substituting  Eq.  (3.4)  into  Eq.  (3.25)  results  in 


Ms)  =  MO)  - 


U1  *  nl  0n  ^1 


(3.26) 


Equation  (3.26)  therefore  corresponds  to  a  Dirichlet  condition 
for  if  on  S1  .  The  only  quantity  on  the  right-hand  side  of  Eq.  (3.26) 
which  is  still  unknown  at  this  time  is  the  constant  value  MO) >  whose 
determination  is  discussed  below. 

3.2.1  Determination  of  the  Arbitrary  Constant 
To  determine  the  constant  MO),  consider  Figure  3.5. 


Fig.  3.5  Object  bounded  by  the  curve  . 

Let  s  *  0  at  an  arbitrarily  chosen  point  P  on  S^,  and  let  Q  be  another 
point  on  S^,  a  distance  s  away  from  P  and  which  is  measured  along  in 
a  counter-clockwise  direction.  Furthermore  let  the  stream  function  be 
denoted  by  ^(0)  and  \p( s)  at  P  and  Q,  respectively. 

In  considering  the  Dirichlet  condition  (3.26),  let  it  be 
assumed  that  ^(0)  at  P  is  given  by  the  constant  a^.  It  is  clear  that, 
through  Eq.  (3.26),  iHs)  is  completely  described  on  provided  the 
constant  a^  is  known.  This  constant  may  be  evaluated  as  follows. 


Consider  Figure  3.6. 


uniform 

stream 


Q(s) :  iKs) 


P(s=0):  i^(0) 


uniform 

stream 


Fig,  3.6  Object  moving  in  an  otherwise  uniform  stream. 


Let  a  new  stream  function  be  defined  by 


^1  +  al^2 


where  ^  and  are  defined  by  the  following  subproblems. 


1.  Find  \p  such  that 


(3.27) 


V2^  =0  in  D  , 


subject  to  the  Dirichlet  condition 


(3.28) 


iK  *  Vx  -  Uy  on  S 


0 


(3.29a) 


r 
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or  the  nonhomogeneous  Neumann  condition 


t —  “  -  (q  x  nA)  .  k  on  SA 

3n  uniform  0  0 


and  the  Dirichlet  condition 


U1  *  nl  °n  ^1 


(3.29b) 


(3.30) 


2.  Find  ^  sucb  that 


V2^  =0  in  D  , 


subject  to  the  Dirichlet  conditions 


^2  =  0  on  SQ 


\f>2  =  1  on  si 


(3.31) 


(3.32) 

(3.33) 


From  Eq.  (3.27),  and  the  subproblems  1.  and  2.  as  defined  above,  this 
new  stream  function  tp*  is  seen  to  be  the  solution  to 


V2ii>*  -  0  in  D 


subject  to  the  Dirichlet  condition 


(3.34) 


^*  =  Vx  -  Uy  on  Sr 


or  the  nonhomogeneous  Neumann  condition 


(3.35a) 


d\j>* 

9n 


(quniform  x  n0)  *  k  on  S0 


(3.35b) 


_  .  J 
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and  the  Dirichlet  condition 


♦*  “  ax 


nx  dS 


on 


(3.36) 


0 

Comparing  Eqs.  (3.34),  (3.35a),  (3.35b),  and  (3.36)  with  Eqs.  (2.12b), 
(3.16),  (3.17b),  and  (3.26),  and  noting  that  tj>(0)  =  a^,  it  is  clear  that 
\p*  as  defined  by  Eq.  (3.27), 


+  ai^2  *  0.27) 

is  the  solution  ip  to  the  problem  at  hand. 

Subproblems  1.  and  2.  are  completely  defined  and  hence  can  be 

solved  using  the  finite  element  method.  The  constant  a^  in  Eq.  (3.27) 

is  determined  as  follows.  Choose  an  arbitrary  point  in  D,  however 

close  to  Sq,  say  B^.  Then,  since  the  flow  is  uniform  at  Sq,  it  is 

expected  that  the  flow  is  still  uniform  at  this  newly  chosen  point  B.  , 

in 

provided  it  is  close  to  S^.  Consequently  the  value  for  the  stream 

function  at  B,  must  be  the  same  as  that  for  the  point  B  _  in  the  outer 
in  r  out 

field  near  Sn,  provided  that  B,  and  B  lie  on  the  same  streamline, 

0  in  out  * 


see  Figure  3.7. 


Fig.  3.7  Streamlines  crossing  the  boundary  S^. 


From  Eq.  (3.27),  and  dropping  the  superscript  *, 
=  •>l'i  +  2^2 


(3. 


and  applying  this  equation  to  the  point  ®^n*  there  obtains 


■  *1<V 


<3. 


Substituting  the  finite  element  solutions  for  and  at  t^e  Point 
Bin  into  Eq.  (3.38),  and  using  the  known  value  of  ^(Bout)  *n 

*<V  •  *<Bout>  -  <3- 


yields  after  rearrangement  of  Eq.  (3.38) 


26. 


♦(Bout>  -  W 

*2<V 


(3.40) 


Equation  (3.40)  yields  a,  provided  that  ty0(B.  )  #  0.  It  is  noted  from 

i  i  in 

Subproblem  2.  that  ^  is  different  from  zero  in  D,  and  hence  ^2^Bin^  ^  0 

for  B,  not  on  SA. 
in  0 

In  passing,  it  is  pointed  out  that  in  the  above  analysis  only 
one  constant  may  be  chosen  arbitrarily,  i.e.,  in  Eq.  (3.16)  C  was  chosen 
to  be  zero.  Consequently  the  need  for  the  evaluation  of  the  constant 
is  apparent.  However,  if  the  nonhomogeneous  Neumann  boundary 
condition  (3.17b)  on  were  chosen  rather  than  the  Dirichlet  condition 
(3.16),  then  there  would  be  no  need  to  perform  the  above  analysis  to 
obtain  a^.  However,  if  more  than  one  object  is  considered  the  choice 
of  Eq.  (3.17)  rather  than  (3.16)  does  not  remove  this  necessity  to 
determine  another  constant,  since  an  additional  constant  defined  by 
b^  *  *(0)  on  $2  must  be  evaluated  in  the  manner  outlined  above. 

It  is  interesting  to  consider  how  the  value  of  »  ^(0)  at 
the  fixed  boundary  point  P  changes  from  its  value  at  time  t  to  its  value 
at  time  t  +  At.  Since  for  this  fixed  point  P  on  the  boundary  of  the 
object  ,  s  was  chosen  to  be  zero,  it  retains  this  value  (s  *  0)  with 
changes  of  time. 

In  general,  from 

i>  -  <p(s,t)  ,  (3.22) 


the  differential  of  (s » t )  is  given  by 


dij)(s,t)  =  -ds  +  -dt 


At  the  point  P,  since  s  =  0  for  all  times,  it  follows  that 


d'l'p(t)  =  d*(0,t)  =  dt 


and  consequently 


^p(t+At)  =  ^(0, t+At)  =  <K0,  t)  + 


3^(0, t) 


Equation  (3.42)  can  be  used  if  the  term 


is  known. 


Alternatively,  i/>p(t+At)  can  be  determined  by  repeating  the  previous 
analysis,  as  follows.  After  the  increment  of  time  At,  the  object  A^ 
occupies  its  new  location,  which  is  completely  determined,  presuming 
u^(t)  is  known.  For  this  new  position  and  this  new  time  t  +  At ,  the 
earlier  analysis  can  be  repeated  with  a  new  constant  a^  being  determined 
as  before.  This  procedure  can  then  be  repeated  for  subsequent  time 
intervals . 


4.  THE  FINITE  ELEMENT  SOLUTION  FOR  UNSTEADY  FLOW 


As  indicated  in  the  previous  sections,  both  a  velocity  potential 
<j>  and  a  stream  function  exist  for  unsteady,  irrotational  flow,  which 
satisfy  Laplace’s  equation  if  the  fluid  is  incompressible. 

The  finite  element  method  can  thus  be  applied  in  the  usual 


manner  [2],  except  that  the  boundary  conditions  will  be  of  the  Dirichlet 
as  well  as  the  nonhomogeneous  Neumann  type  when  the  flow  field  involves 
moving  boundaries,  see  Section  3.  This  means  that  the  mathematical 


statement  of  the  problem  for  either  the  velocity  potential  4>  or  the 
stream  function  at  any  given  instant  of  time  reduces  to  the  follow 
ing  (see  Figure  4.1). 


h. 


X 


0  x 

Fig.  4.1  Domain  D  enclosed  by  the  surfaces  Sq  and  S^. 

Obtain  that  function  uCx.y.t^),  for  t^  arbitrary,  which 
satisfies  Laplace’s  equation  in  the  domain  D,  .e.. 


V^Cxjy ,tQ)  =  0  in  D 


subject  to  the  following  boundary  conditions: 
1.  u(x,y,tQ)  -  g(x,y,tQ)  on  Sp 


or 

3u(x,y,tQ) 


on 


S 


0 


(4. 


(4. 


3n 


h(x,y,tft) 


(4. 
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3u(x,y,t  ) 

2.  - ^ -  -  z(x,y,tQ)  on  Sx  ,  (4.3) 

where  g(x,y,tg)  and  h(x,y,tQ)  are  prescribed  functions  of  position 
(and  time^)  along  the  fixed  (i.e.,  with  respect  to  the  fixed 
reference  frame  Oxy)  boundary  S^,  z(x,y,tg)  is  a  prescribed  function 
of  position  and  time  t^  along  S^.  Although  the  boundary  moves 

with  time,  its  location  may  be  considered  fixed  and  completely  pre¬ 
scribed  with  respect  to  the  fixed  reference  frame  Oxy  at  any 
instant  of  time  and  hence  also  at  t^. 

The  problem  stated  above  may  be  reformulated  in  terms  of  the 
calculus  of  variations  as  follows.  Obtain  that  function  v(x,y,t^),  which 
minimizes  the  functional 

X(v)  =  f  M(|J)2  +  (|p02]dD  -  A  |  h(x,y,t0)vdS0  -  j  z(x,y,t())vdS1 

»  S0  S1 

(4.4) 

where  A  is  a  constant,  either  0  or  1  in  value,  depending  on  which  form¬ 
ulation  is  chosen,  i.e.,  A  -  0  for  the  problem  defined  by  Eqs.  (4,1), 

(4.2a),  (4.3),  and  A  =  1  if  the  governing  equations  (4,1),  (4,2b),  and 
(4.3)  are  used.  In  the  functional  (4.4),  the  function  v  =  v(x,y,t^)  must 
belong  to  the  class  of  admissible  functions,  i.e.,  it  must  be  continuous 
and  have  piecewise  continuous  first  derivatives  in  D,  and  furthermore 
it  must  satisfy  the  principal  or  Dirichlet  boundary  conditions.  This 
latter  requirement,  of  course,  only  applies  if  the  governing  equations 

t  If  the  case  of  varying  free-stream  velocity  is  included. 
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(4.1),  (4.2a),  and  (4.3)  are  considered. 

It  can  be  shown  [2]  that  that  function  v(x,y,tg)  which 
minimizes  the  given  functional  (4.4),  with  the  proper  choice  of  the 
constant  A,  is  also  a  solution  to  the  field  problem  defined  by 
Eqs.  (4.1),  (4.2a)  or  (4.2b),  and  (4.3).  The  finite  element  method 
as  applied  to  the  functional  (4.4)  has  been  described  in  the  literature 
[2,8,9],  and  will  not  be  considered  any  further  in  this  report. 

4.1  The  Velocity  Potential 

Comparison  of  Eqs.  (2.8b),  (3.8),  (3.13),  and  (3.4),  with 
Eqs.  (4.1),  (4.2a),  (4.2b),  and  (4.3)  indicate  that  the  solution  to  the 
velocity  potential  problem  considered  in  Sections  (2.1)  and  (3.1)  can 
be  obtained  using  the  above  analysis  and  making  the  following  choice 
for  the  variables  g,  h,  and  z,  respectively, 

g(x,y,tfl)  *  -  Ux  -  Vy  on  SQ  ,  (4.5a] 

h(x,y,to)  -  -  5unifom  •  50  on  S0  ,  (4.5b; 

zU.y.tg)  -  -  ^(Cq)  .  nx  on  .  (4.5c] 

Substituting  the  expressions  (4.5b),  (4.5c)  into  the  functional  (4.4), 
and  applying  the  finite  element  technique  to  this  functional  yields  the 
solution  for  <j>  ■  <J>(x,y,tn)  . 
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4.2  The  Stream  Function 

Comparison  of  Eqs,  (2.12b),  (3.16),  (3.17b),  and  (3.26),  with 
Eqs.  (4.1),  (4.2a),  (4.2b),  and  (4.3)  indicate  that  the  solution  to 
the  stream  function  problem  considered  in  Sections  (2.2)  and  (3.2)  can 
be  obtained  using  the  above  analysis  and  making  the  following  choice 


for  the  variables  g, 

g(x,y,tQ) 

h,  and  z,  respectively, 

=  Vx  -  Uy  on  , 

(4.6a) 

h(x,y,tQ) 

*  -  (q  x  nrt)  .  k  on  S. 

^uniform  0  0 

I 

(4.6b) 

z(x,y,t0) 

=  a  -  1 

S 

* 

u^  .  n^  ds  on  S^ 

♦ 

(4.6c) 

j 

0 

Substituting  the  expressions  (4.6b)  and  (4.6c)  into  the  functional  (4.4), 
and  applying  the  finite  element  technique  to  this  functional  yields  the 
solution  for  ip  *  ^(x,y,t^), 

5.  ENTRAINED  MASS  OR  INERTIA 

In  moving  body  problems  such  as  described  in  the  previous 
sections,  the  pressure  can  be  calculated  from  Eq.  (2.9), 

£  -  f£  +  \  (V<f>)2  =  f(t)  ,  (2.9) 

which  can  be  rearranged  as 

p  =  -  y  (V<J>)2  +  p  +  pf(t)  ,  (5.1) 

or  as 

=  p  +  p  +  P« 

*v  ra  r0 


P 


(5.2) 


where  the  first  right-hand-side  term  is  interpreted  as  the  pressure 
related  to  velocity,  and  the  second  term  as  the  pressure  related  to 
acceleration  since  it  disappears  in  the  absence  of  acceleration. 

The  integral  of  p^  around  the  body  surface  gives  the  force 
corresponding  to  steady-state  conditions,  whereas  the  integral  of  p 

a 

gives  the  'entrained  mass'  force.  For  a  non-lifting  body  in  an 
irrotational  stream,  the  latter  acceleration  force  divided  by  the  in¬ 
stantaneous  acceleration  gives  a  constant  dependent  on  the  shape  and 
orientation  of  the  body  known  as  the  entrained  mass. 

The  finite  element  approaches  outlined  earlier  in  this  report 
can  be  used  to  solve  entrained  mass  problems  since  Eq.  (5.1)  can  be 
evaluated  at  a  given  point  on  the  body  surface  at  a  particular  time  t^ 
with  V<£  and  being  obtained  as  follows: 

(1)  7<J)  can  be  obtained  from  the  finite  element  solution  at  the  given 
time  t-^  by  either 

a)  using  as  the  unknown  and  deriving  from  the  resultant 
solution,  or 

b)  using  a  higher-order  shape  function  solution  with  4>  and  the 
elements  of  V<p  as  unknowns, 

(2)  *!£  at  t^  can  be  derived  by  using  the  solutions  for  t  ®  t^  -  At 
and  t  *  t^  and  extrapolating  from  at  (t^At)  using  a  suitable 
formula. 

Integration  of  the  pressure  around  the  surface  allows  the  acceleration 
force  and  hence  the  entrained  mass  or  inertia  to  be  derived. 

The  stream  function  solutions  given  previously  can  be  used 
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